
R version 3.1.0 (2014-04-10) -- "Spring Dance"
Copyright (C) 2014 The R Foundation for Statistical Computing
Platform: x86_64-unknown-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

  Natural language support but running in an English locale

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> set.seed(51)
> 
> load(file="~/PM10/pm10_withps1990.Rda")
> ######## Prune observations that are out of range ###########
> dat=subset(dat, outofrange==0)
> dat$pscat = as.factor(dat$pscat)
> 
> ######## Set the Working Directory so that output is saved in the right place ######## 
> wd <- "~/PM10/20150502_big_1990/"
> setwd(wd)
> 
> ##############################################################
> ################ Specify the Array Parameters ################
> ##############################################################
> 
> outcomenames = c("Total_death.2001", "ALL_CVD.2001", "Respiratory.2001")  ### The outcomes that will be considered
> denominators = c("Tot_Beneficiary.2001", "pyear.2001", "pyear.2001")
> outfiles = c("mort", "cvd", "resp")                                       ### To append the output R object name
> 
> 
> ##############################################################
> ###################### Array Starts Here #####################
> ##############################################################
> j = as.numeric(Sys.getenv('SLURM_ARRAY_TASK_ID'))
> print(j)
[1] 1
> 
> ######## Choose which health outcome to use ######## 
> dat$h = dat[, outcomenames[j]]                                             ### Set h (outcome variable)
> dat$denom = dat[, denominators[j]]                                         ### Set denom (denominator variable)
> 
> 
> ######## Choose how to transform pollution ######## 
> dat$y = log(dat$pm10fu)    #dat$fu - dat$base
> dat$ytemp=dat$y
> dat$ytemp[is.na(dat$ytemp)]= mean(dat$y, na.rm=TRUE)
> 
> ######## Choose how to specify the models ######## 
> ### Big Models ###
> mh0=glm(h~pscat + pm10base1990 + Median_income + HS_rate + Urban_rate + Migration_5_year_rate + Current_Smoking_rate_weighted + mean_age.2001 + Female_rate.2001 + White_rate.2001 + Black_rate.2001 + housedens + temperature + logpop +  ytemp + offset(log(denom)), 
+         family=poisson(link="log"), data=subset(dat, a==0))
> mh1=glm(h~pscat + pm10base1990 + Median_income + HS_rate + Urban_rate + Migration_5_year_rate + Current_Smoking_rate_weighted + mean_age.2001 + Female_rate.2001 + White_rate.2001 + Black_rate.2001 + housedens + temperature + logpop +  ytemp + offset(log(denom)), 
+         family=poisson(link="log"), data=subset(dat, a==1))
> 
> 
> library(HEIfunctions)
> 
> q=2
> n=dim(dat)[[1]]
> nburn= 0 #5000 the function only outputs iters > nburn
> nsamp= 50000 ## nburnin + niter
> thin=1
> alpha0prop=summary(mh0)$cov.unscaled*1.1
> alpha1prop=summary(mh1)$cov.unscaled*.9
> ypropsd=rep(.75, n*q)
> 
> #### Proposals for log(pm) ######
> ypropsd[seq(1,n*q,q)] = log(dat$pm10base1990_1994)*.3  
> ypropsd[seq(2,n*q,q)] = log(dat$pm10base1990_1994)*.3 
> 
> #### Proposals for pmfu - pmbase ######
> # ypropsd[seq(1,n*q,q)] = dat$pm10base1990_1994*.5
> # ypropsd[seq(2,n*q,q)] = dat$pm10base1990_1994*.5 
> 
> coords   	<- cbind(dat$Longitude, dat$Latitude)
> phistart 	<- makephi(coords, 10)
> philower 	<- makephi(coords, 4) 
> phiupper 	<- makephi(coords, 30)
> 
> tuning = list(A = 0.1, psi = 0.2, theta = 1, alpha0=.09*alpha0prop, alpha1=.09*alpha1prop, Y=ypropsd)
> prior = list(KIG = rep(.5,2), psi = rep(.5,2), theta1 = rep(philower,2), theta2 = rep(phiupper,2))
> starting = list(B=NULL, A = NULL, psi = NULL, theta = phistart, alpha0=mh0$coef, alpha1=mh1$coef)
> 
> 
> pollutionformula = as.formula( paste("y ~", paste(all.vars(mh0$formula)[!(all.vars(mh0$formula) %in% c("h", "ytemp", "denom"))], collapse="+"))  )
> healthformula = as.formula( paste("h ~", paste(all.vars(mh0$formula)[!(all.vars(mh0$formula) %in% c("h", "ytemp", "denom"))], collapse="+"))  )
> 
> starttime=proc.time()
> pstratamod=principalstratmod(
+ formula = pollutionformula, data=dat, trt="a", coords, 
+ formula_h = healthformula, denom="denom", nsamp, thin, tuning, prior, starting,
+ outputbinsize = 500, outputfilename = paste("pstratamod_temp_", outfiles[j], ".Rda", sep=""))
Loading required package: spBayes
Loading required package: coda
Loading required package: lattice
Loading required package: magic
Loading required package: abind
Loading required package: Formula
Loading required package: corpcor
Loading required package: MASS
[1] "Note: function does not accept missing data in covariates."
[1] "Iteration: 200"
 [1] "Accepted:" "0.27"      "0.415"     "0.27"      "0.415"     "0.27"     
 [7] "0.415"     "0.185"     "0.185"     "0.185"     "0.185"     "0.185"    
[13] "0.185"     "0.185"     "0.185"     "0.185"     "0.185"     "0.185"    
[19] "0.185"     "0.185"     "0.185"     "0.185"     "0.185"     "0.185"    
[25] "0.185"     "0.185"     "0.345"     "0.345"     "0.345"     "0.345"    
[31] "0.345"     "0.345"     "0.345"     "0.345"     "0.345"     "0.345"    
[37] "0.345"     "0.345"     "0.345"     "0.345"     "0.345"     "0.345"    
[43] "0.345"     "0.345"     "0.345"     "0.3"      
[1] "Number of not PD:" "0"                
[1] "Iteration: 400"
 [1] "Accepted:" "0.295"     "0.38"      "0.295"     "0.38"      "0.295"    
 [7] "0.38"      "0.198"     "0.198"     "0.198"     "0.198"     "0.198"    
[13] "0.198"     "0.198"     "0.198"     "0.198"     "0.198"     "0.198"    
[19] "0.198"     "0.198"     "0.198"     "0.198"     "0.198"     "0.198"    
[25] "0.198"     "0.198"     "0.33"      "0.33"      "0.33"      "0.33"     
[31] "0.33"      "0.33"      "0.33"      "0.33"      "0.33"      "0.33"     
[37] "0.33"      "0.33"      "0.33"      "0.33"      "0.33"      "0.33"     
[43] "0.33"      "0.33"      "0.33"      "0.301"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 600"
 [1] "Accepted:" "0.278"     "0.368"     "0.278"     "0.368"     "0.278"    
 [7] "0.368"     "0.187"     "0.187"     "0.187"     "0.187"     "0.187"    
[13] "0.187"     "0.187"     "0.187"     "0.187"     "0.187"     "0.187"    
[19] "0.187"     "0.187"     "0.187"     "0.187"     "0.187"     "0.187"    
[25] "0.187"     "0.187"     "0.305"     "0.305"     "0.305"     "0.305"    
[31] "0.305"     "0.305"     "0.305"     "0.305"     "0.305"     "0.305"    
[37] "0.305"     "0.305"     "0.305"     "0.305"     "0.305"     "0.305"    
[43] "0.305"     "0.305"     "0.305"     "0.297"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 800"
 [1] "Accepted:" "0.271"     "0.352"     "0.271"     "0.352"     "0.271"    
 [7] "0.352"     "0.196"     "0.196"     "0.196"     "0.196"     "0.196"    
[13] "0.196"     "0.196"     "0.196"     "0.196"     "0.196"     "0.196"    
[19] "0.196"     "0.196"     "0.196"     "0.196"     "0.196"     "0.196"    
[25] "0.196"     "0.196"     "0.295"     "0.295"     "0.295"     "0.295"    
[31] "0.295"     "0.295"     "0.295"     "0.295"     "0.295"     "0.295"    
[37] "0.295"     "0.295"     "0.295"     "0.295"     "0.295"     "0.295"    
[43] "0.295"     "0.295"     "0.295"     "0.295"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 1000"
 [1] "Accepted:" "0.274"     "0.339"     "0.274"     "0.339"     "0.274"    
 [7] "0.339"     "0.204"     "0.204"     "0.204"     "0.204"     "0.204"    
[13] "0.204"     "0.204"     "0.204"     "0.204"     "0.204"     "0.204"    
[19] "0.204"     "0.204"     "0.204"     "0.204"     "0.204"     "0.204"    
[25] "0.204"     "0.204"     "0.294"     "0.294"     "0.294"     "0.294"    
[31] "0.294"     "0.294"     "0.294"     "0.294"     "0.294"     "0.294"    
[37] "0.294"     "0.294"     "0.294"     "0.294"     "0.294"     "0.294"    
[43] "0.294"     "0.294"     "0.294"     "0.293"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 1200"
 [1] "Accepted:" "0.269"     "0.344"     "0.269"     "0.344"     "0.269"    
 [7] "0.344"     "0.196"     "0.196"     "0.196"     "0.196"     "0.196"    
[13] "0.196"     "0.196"     "0.196"     "0.196"     "0.196"     "0.196"    
[19] "0.196"     "0.196"     "0.196"     "0.196"     "0.196"     "0.196"    
[25] "0.196"     "0.196"     "0.289"     "0.289"     "0.289"     "0.289"    
[31] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[37] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[43] "0.289"     "0.289"     "0.289"     "0.292"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 1400"
 [1] "Accepted:" "0.269"     "0.353"     "0.269"     "0.353"     "0.269"    
 [7] "0.353"     "0.192"     "0.192"     "0.192"     "0.192"     "0.192"    
[13] "0.192"     "0.192"     "0.192"     "0.192"     "0.192"     "0.192"    
[19] "0.192"     "0.192"     "0.192"     "0.192"     "0.192"     "0.192"    
[25] "0.192"     "0.192"     "0.289"     "0.289"     "0.289"     "0.289"    
[31] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[37] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[43] "0.289"     "0.289"     "0.289"     "0.29"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 1600"
 [1] "Accepted:" "0.273"     "0.351"     "0.273"     "0.351"     "0.273"    
 [7] "0.351"     "0.192"     "0.192"     "0.192"     "0.192"     "0.192"    
[13] "0.192"     "0.192"     "0.192"     "0.192"     "0.192"     "0.192"    
[19] "0.192"     "0.192"     "0.192"     "0.192"     "0.192"     "0.192"    
[25] "0.192"     "0.192"     "0.293"     "0.293"     "0.293"     "0.293"    
[31] "0.293"     "0.293"     "0.293"     "0.293"     "0.293"     "0.293"    
[37] "0.293"     "0.293"     "0.293"     "0.293"     "0.293"     "0.293"    
[43] "0.293"     "0.293"     "0.293"     "0.289"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 1800"
 [1] "Accepted:" "0.276"     "0.346"     "0.276"     "0.346"     "0.276"    
 [7] "0.346"     "0.196"     "0.196"     "0.196"     "0.196"     "0.196"    
[13] "0.196"     "0.196"     "0.196"     "0.196"     "0.196"     "0.196"    
[19] "0.196"     "0.196"     "0.196"     "0.196"     "0.196"     "0.196"    
[25] "0.196"     "0.196"     "0.293"     "0.293"     "0.293"     "0.293"    
[31] "0.293"     "0.293"     "0.293"     "0.293"     "0.293"     "0.293"    
[37] "0.293"     "0.293"     "0.293"     "0.293"     "0.293"     "0.293"    
[43] "0.293"     "0.293"     "0.293"     "0.288"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 2000"
 [1] "Accepted:" "0.282"     "0.337"     "0.282"     "0.337"     "0.282"    
 [7] "0.337"     "0.198"     "0.198"     "0.198"     "0.198"     "0.198"    
[13] "0.198"     "0.198"     "0.198"     "0.198"     "0.198"     "0.198"    
[19] "0.198"     "0.198"     "0.198"     "0.198"     "0.198"     "0.198"    
[25] "0.198"     "0.198"     "0.288"     "0.288"     "0.288"     "0.288"    
[31] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[37] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[43] "0.288"     "0.288"     "0.288"     "0.288"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 2200"
 [1] "Accepted:" "0.285"     "0.336"     "0.285"     "0.336"     "0.285"    
 [7] "0.336"     "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[13] "0.2"       "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[19] "0.2"       "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[25] "0.2"       "0.2"       "0.287"     "0.287"     "0.287"     "0.287"    
[31] "0.287"     "0.287"     "0.287"     "0.287"     "0.287"     "0.287"    
[37] "0.287"     "0.287"     "0.287"     "0.287"     "0.287"     "0.287"    
[43] "0.287"     "0.287"     "0.287"     "0.287"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 2400"
 [1] "Accepted:" "0.288"     "0.336"     "0.288"     "0.336"     "0.288"    
 [7] "0.336"     "0.198"     "0.198"     "0.198"     "0.198"     "0.198"    
[13] "0.198"     "0.198"     "0.198"     "0.198"     "0.198"     "0.198"    
[19] "0.198"     "0.198"     "0.198"     "0.198"     "0.198"     "0.198"    
[25] "0.198"     "0.198"     "0.288"     "0.288"     "0.288"     "0.288"    
[31] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[37] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[43] "0.288"     "0.288"     "0.288"     "0.287"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 2600"
 [1] "Accepted:" "0.288"     "0.331"     "0.288"     "0.331"     "0.288"    
 [7] "0.331"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[13] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[19] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[25] "0.201"     "0.201"     "0.289"     "0.289"     "0.289"     "0.289"    
[31] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[37] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[43] "0.289"     "0.289"     "0.289"     "0.287"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 2800"
 [1] "Accepted:" "0.29"      "0.331"     "0.29"      "0.331"     "0.29"     
 [7] "0.331"     "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[13] "0.2"       "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[19] "0.2"       "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[25] "0.2"       "0.2"       "0.29"      "0.29"      "0.29"      "0.29"     
[31] "0.29"      "0.29"      "0.29"      "0.29"      "0.29"      "0.29"     
[37] "0.29"      "0.29"      "0.29"      "0.29"      "0.29"      "0.29"     
[43] "0.29"      "0.29"      "0.29"      "0.286"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 3000"
 [1] "Accepted:" "0.289"     "0.332"     "0.289"     "0.332"     "0.289"    
 [7] "0.332"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[13] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[19] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[25] "0.201"     "0.201"     "0.289"     "0.289"     "0.289"     "0.289"    
[31] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[37] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[43] "0.289"     "0.289"     "0.289"     "0.287"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 3200"
 [1] "Accepted:" "0.29"      "0.332"     "0.29"      "0.332"     "0.29"     
 [7] "0.332"     "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[13] "0.2"       "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[19] "0.2"       "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[25] "0.2"       "0.2"       "0.288"     "0.288"     "0.288"     "0.288"    
[31] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[37] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[43] "0.288"     "0.288"     "0.288"     "0.287"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 3400"
 [1] "Accepted:" "0.288"     "0.333"     "0.288"     "0.333"     "0.288"    
 [7] "0.333"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[13] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[19] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[25] "0.201"     "0.201"     "0.288"     "0.288"     "0.288"     "0.288"    
[31] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[37] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[43] "0.288"     "0.288"     "0.288"     "0.286"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 3600"
 [1] "Accepted:" "0.288"     "0.332"     "0.288"     "0.332"     "0.288"    
 [7] "0.332"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[13] "0.202"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[19] "0.202"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[25] "0.202"     "0.202"     "0.289"     "0.289"     "0.289"     "0.289"    
[31] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[37] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[43] "0.289"     "0.289"     "0.289"     "0.286"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 3800"
 [1] "Accepted:" "0.289"     "0.332"     "0.289"     "0.332"     "0.289"    
 [7] "0.332"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[13] "0.203"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[19] "0.203"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[25] "0.203"     "0.203"     "0.289"     "0.289"     "0.289"     "0.289"    
[31] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[37] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[43] "0.289"     "0.289"     "0.289"     "0.286"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 4000"
 [1] "Accepted:" "0.288"     "0.329"     "0.288"     "0.329"     "0.288"    
 [7] "0.329"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[13] "0.202"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[19] "0.202"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[25] "0.202"     "0.202"     "0.291"     "0.291"     "0.291"     "0.291"    
[31] "0.291"     "0.291"     "0.291"     "0.291"     "0.291"     "0.291"    
[37] "0.291"     "0.291"     "0.291"     "0.291"     "0.291"     "0.291"    
[43] "0.291"     "0.291"     "0.291"     "0.285"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 4200"
 [1] "Accepted:" "0.289"     "0.327"     "0.289"     "0.327"     "0.289"    
 [7] "0.327"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[13] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[19] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[25] "0.201"     "0.201"     "0.291"     "0.291"     "0.291"     "0.291"    
[31] "0.291"     "0.291"     "0.291"     "0.291"     "0.291"     "0.291"    
[37] "0.291"     "0.291"     "0.291"     "0.291"     "0.291"     "0.291"    
[43] "0.291"     "0.291"     "0.291"     "0.285"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 4400"
 [1] "Accepted:" "0.289"     "0.326"     "0.289"     "0.326"     "0.289"    
 [7] "0.326"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[13] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[19] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[25] "0.201"     "0.201"     "0.291"     "0.291"     "0.291"     "0.291"    
[31] "0.291"     "0.291"     "0.291"     "0.291"     "0.291"     "0.291"    
[37] "0.291"     "0.291"     "0.291"     "0.291"     "0.291"     "0.291"    
[43] "0.291"     "0.291"     "0.291"     "0.284"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 4600"
 [1] "Accepted:" "0.29"      "0.324"     "0.29"      "0.324"     "0.29"     
 [7] "0.324"     "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[13] "0.2"       "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[19] "0.2"       "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[25] "0.2"       "0.2"       "0.293"     "0.293"     "0.293"     "0.293"    
[31] "0.293"     "0.293"     "0.293"     "0.293"     "0.293"     "0.293"    
[37] "0.293"     "0.293"     "0.293"     "0.293"     "0.293"     "0.293"    
[43] "0.293"     "0.293"     "0.293"     "0.283"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 4800"
 [1] "Accepted:" "0.294"     "0.323"     "0.294"     "0.323"     "0.294"    
 [7] "0.323"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[13] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[19] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[25] "0.201"     "0.201"     "0.291"     "0.291"     "0.291"     "0.291"    
[31] "0.291"     "0.291"     "0.291"     "0.291"     "0.291"     "0.291"    
[37] "0.291"     "0.291"     "0.291"     "0.291"     "0.291"     "0.291"    
[43] "0.291"     "0.291"     "0.291"     "0.283"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 5000"
 [1] "Accepted:" "0.292"     "0.325"     "0.292"     "0.325"     "0.292"    
 [7] "0.325"     "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[13] "0.2"       "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[19] "0.2"       "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[25] "0.2"       "0.2"       "0.289"     "0.289"     "0.289"     "0.289"    
[31] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[37] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[43] "0.289"     "0.289"     "0.289"     "0.283"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 5200"
 [1] "Accepted:" "0.294"     "0.327"     "0.294"     "0.327"     "0.294"    
 [7] "0.327"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[13] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[19] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[25] "0.201"     "0.201"     "0.289"     "0.289"     "0.289"     "0.289"    
[31] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[37] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[43] "0.289"     "0.289"     "0.289"     "0.283"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 5400"
 [1] "Accepted:" "0.297"     "0.326"     "0.297"     "0.326"     "0.297"    
 [7] "0.326"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[13] "0.202"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[19] "0.202"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[25] "0.202"     "0.202"     "0.291"     "0.291"     "0.291"     "0.291"    
[31] "0.291"     "0.291"     "0.291"     "0.291"     "0.291"     "0.291"    
[37] "0.291"     "0.291"     "0.291"     "0.291"     "0.291"     "0.291"    
[43] "0.291"     "0.291"     "0.291"     "0.283"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 5600"
 [1] "Accepted:" "0.297"     "0.325"     "0.297"     "0.325"     "0.297"    
 [7] "0.325"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[13] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[19] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[25] "0.201"     "0.201"     "0.291"     "0.291"     "0.291"     "0.291"    
[31] "0.291"     "0.291"     "0.291"     "0.291"     "0.291"     "0.291"    
[37] "0.291"     "0.291"     "0.291"     "0.291"     "0.291"     "0.291"    
[43] "0.291"     "0.291"     "0.291"     "0.283"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 5800"
 [1] "Accepted:" "0.295"     "0.324"     "0.295"     "0.324"     "0.295"    
 [7] "0.324"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[13] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[19] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[25] "0.201"     "0.201"     "0.291"     "0.291"     "0.291"     "0.291"    
[31] "0.291"     "0.291"     "0.291"     "0.291"     "0.291"     "0.291"    
[37] "0.291"     "0.291"     "0.291"     "0.291"     "0.291"     "0.291"    
[43] "0.291"     "0.291"     "0.291"     "0.283"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 6000"
 [1] "Accepted:" "0.295"     "0.326"     "0.295"     "0.326"     "0.295"    
 [7] "0.326"     "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[13] "0.2"       "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[19] "0.2"       "0.2"       "0.2"       "0.2"       "0.2"       "0.2"      
[25] "0.2"       "0.2"       "0.292"     "0.292"     "0.292"     "0.292"    
[31] "0.292"     "0.292"     "0.292"     "0.292"     "0.292"     "0.292"    
[37] "0.292"     "0.292"     "0.292"     "0.292"     "0.292"     "0.292"    
[43] "0.292"     "0.292"     "0.292"     "0.283"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 6200"
 [1] "Accepted:" "0.297"     "0.324"     "0.297"     "0.324"     "0.297"    
 [7] "0.324"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[13] "0.202"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[19] "0.202"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[25] "0.202"     "0.202"     "0.289"     "0.289"     "0.289"     "0.289"    
[31] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[37] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[43] "0.289"     "0.289"     "0.289"     "0.283"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 6400"
 [1] "Accepted:" "0.297"     "0.324"     "0.297"     "0.324"     "0.297"    
 [7] "0.324"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[13] "0.202"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[19] "0.202"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[25] "0.202"     "0.202"     "0.29"      "0.29"      "0.29"      "0.29"     
[31] "0.29"      "0.29"      "0.29"      "0.29"      "0.29"      "0.29"     
[37] "0.29"      "0.29"      "0.29"      "0.29"      "0.29"      "0.29"     
[43] "0.29"      "0.29"      "0.29"      "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 6600"
 [1] "Accepted:" "0.298"     "0.324"     "0.298"     "0.324"     "0.298"    
 [7] "0.324"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[13] "0.202"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[19] "0.202"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[25] "0.202"     "0.202"     "0.29"      "0.29"      "0.29"      "0.29"     
[31] "0.29"      "0.29"      "0.29"      "0.29"      "0.29"      "0.29"     
[37] "0.29"      "0.29"      "0.29"      "0.29"      "0.29"      "0.29"     
[43] "0.29"      "0.29"      "0.29"      "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 6800"
 [1] "Accepted:" "0.298"     "0.323"     "0.298"     "0.323"     "0.298"    
 [7] "0.323"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[13] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[19] "0.201"     "0.201"     "0.201"     "0.201"     "0.201"     "0.201"    
[25] "0.201"     "0.201"     "0.289"     "0.289"     "0.289"     "0.289"    
[31] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[37] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[43] "0.289"     "0.289"     "0.289"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 7000"
 [1] "Accepted:" "0.297"     "0.322"     "0.297"     "0.322"     "0.297"    
 [7] "0.322"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[13] "0.203"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[19] "0.203"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[25] "0.203"     "0.203"     "0.286"     "0.286"     "0.286"     "0.286"    
[31] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[37] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[43] "0.286"     "0.286"     "0.286"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 7200"
 [1] "Accepted:" "0.297"     "0.322"     "0.297"     "0.322"     "0.297"    
 [7] "0.322"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[13] "0.203"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[19] "0.203"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[25] "0.203"     "0.203"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 7400"
 [1] "Accepted:" "0.298"     "0.321"     "0.298"     "0.321"     "0.298"    
 [7] "0.321"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[13] "0.203"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[19] "0.203"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[25] "0.203"     "0.203"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 7600"
 [1] "Accepted:" "0.297"     "0.321"     "0.297"     "0.321"     "0.297"    
 [7] "0.321"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[13] "0.202"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[19] "0.202"     "0.202"     "0.202"     "0.202"     "0.202"     "0.202"    
[25] "0.202"     "0.202"     "0.286"     "0.286"     "0.286"     "0.286"    
[31] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[37] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[43] "0.286"     "0.286"     "0.286"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 7800"
 [1] "Accepted:" "0.298"     "0.32"      "0.298"     "0.32"      "0.298"    
 [7] "0.32"      "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[13] "0.203"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[19] "0.203"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[25] "0.203"     "0.203"     "0.286"     "0.286"     "0.286"     "0.286"    
[31] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[37] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[43] "0.286"     "0.286"     "0.286"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 8000"
 [1] "Accepted:" "0.298"     "0.32"      "0.298"     "0.32"      "0.298"    
 [7] "0.32"      "0.204"     "0.204"     "0.204"     "0.204"     "0.204"    
[13] "0.204"     "0.204"     "0.204"     "0.204"     "0.204"     "0.204"    
[19] "0.204"     "0.204"     "0.204"     "0.204"     "0.204"     "0.204"    
[25] "0.204"     "0.204"     "0.286"     "0.286"     "0.286"     "0.286"    
[31] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[37] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[43] "0.286"     "0.286"     "0.286"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 8200"
 [1] "Accepted:" "0.298"     "0.32"      "0.298"     "0.32"      "0.298"    
 [7] "0.32"      "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[13] "0.203"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[19] "0.203"     "0.203"     "0.203"     "0.203"     "0.203"     "0.203"    
[25] "0.203"     "0.203"     "0.288"     "0.288"     "0.288"     "0.288"    
[31] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[37] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[43] "0.288"     "0.288"     "0.288"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 8400"
 [1] "Accepted:" "0.298"     "0.32"      "0.298"     "0.32"      "0.298"    
 [7] "0.32"      "0.204"     "0.204"     "0.204"     "0.204"     "0.204"    
[13] "0.204"     "0.204"     "0.204"     "0.204"     "0.204"     "0.204"    
[19] "0.204"     "0.204"     "0.204"     "0.204"     "0.204"     "0.204"    
[25] "0.204"     "0.204"     "0.288"     "0.288"     "0.288"     "0.288"    
[31] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[37] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[43] "0.288"     "0.288"     "0.288"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 8600"
 [1] "Accepted:" "0.298"     "0.32"      "0.298"     "0.32"      "0.298"    
 [7] "0.32"      "0.204"     "0.204"     "0.204"     "0.204"     "0.204"    
[13] "0.204"     "0.204"     "0.204"     "0.204"     "0.204"     "0.204"    
[19] "0.204"     "0.204"     "0.204"     "0.204"     "0.204"     "0.204"    
[25] "0.204"     "0.204"     "0.288"     "0.288"     "0.288"     "0.288"    
[31] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[37] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[43] "0.288"     "0.288"     "0.288"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 8800"
 [1] "Accepted:" "0.299"     "0.319"     "0.299"     "0.319"     "0.299"    
 [7] "0.319"     "0.204"     "0.204"     "0.204"     "0.204"     "0.204"    
[13] "0.204"     "0.204"     "0.204"     "0.204"     "0.204"     "0.204"    
[19] "0.204"     "0.204"     "0.204"     "0.204"     "0.204"     "0.204"    
[25] "0.204"     "0.204"     "0.288"     "0.288"     "0.288"     "0.288"    
[31] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[37] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[43] "0.288"     "0.288"     "0.288"     "0.283"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 9000"
 [1] "Accepted:" "0.299"     "0.318"     "0.299"     "0.318"     "0.299"    
 [7] "0.318"     "0.205"     "0.205"     "0.205"     "0.205"     "0.205"    
[13] "0.205"     "0.205"     "0.205"     "0.205"     "0.205"     "0.205"    
[19] "0.205"     "0.205"     "0.205"     "0.205"     "0.205"     "0.205"    
[25] "0.205"     "0.205"     "0.288"     "0.288"     "0.288"     "0.288"    
[31] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[37] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[43] "0.288"     "0.288"     "0.288"     "0.283"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 9200"
 [1] "Accepted:" "0.298"     "0.318"     "0.298"     "0.318"     "0.298"    
 [7] "0.318"     "0.206"     "0.206"     "0.206"     "0.206"     "0.206"    
[13] "0.206"     "0.206"     "0.206"     "0.206"     "0.206"     "0.206"    
[19] "0.206"     "0.206"     "0.206"     "0.206"     "0.206"     "0.206"    
[25] "0.206"     "0.206"     "0.288"     "0.288"     "0.288"     "0.288"    
[31] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[37] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[43] "0.288"     "0.288"     "0.288"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 9400"
 [1] "Accepted:" "0.299"     "0.317"     "0.299"     "0.317"     "0.299"    
 [7] "0.317"     "0.207"     "0.207"     "0.207"     "0.207"     "0.207"    
[13] "0.207"     "0.207"     "0.207"     "0.207"     "0.207"     "0.207"    
[19] "0.207"     "0.207"     "0.207"     "0.207"     "0.207"     "0.207"    
[25] "0.207"     "0.207"     "0.288"     "0.288"     "0.288"     "0.288"    
[31] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[37] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[43] "0.288"     "0.288"     "0.288"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 9600"
 [1] "Accepted:" "0.3"       "0.319"     "0.3"       "0.319"     "0.3"      
 [7] "0.319"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[13] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[19] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[25] "0.208"     "0.208"     "0.289"     "0.289"     "0.289"     "0.289"    
[31] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[37] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[43] "0.289"     "0.289"     "0.289"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 9800"
 [1] "Accepted:" "0.3"       "0.318"     "0.3"       "0.318"     "0.3"      
 [7] "0.318"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.288"     "0.288"     "0.288"     "0.288"    
[31] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[37] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[43] "0.288"     "0.288"     "0.288"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 10000"
 [1] "Accepted:" "0.3"       "0.318"     "0.3"       "0.318"     "0.3"      
 [7] "0.318"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[13] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[19] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[25] "0.208"     "0.208"     "0.289"     "0.289"     "0.289"     "0.289"    
[31] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[37] "0.289"     "0.289"     "0.289"     "0.289"     "0.289"     "0.289"    
[43] "0.289"     "0.289"     "0.289"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 10200"
 [1] "Accepted:" "0.301"     "0.318"     "0.301"     "0.318"     "0.301"    
 [7] "0.318"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[13] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[19] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[25] "0.208"     "0.208"     "0.288"     "0.288"     "0.288"     "0.288"    
[31] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[37] "0.288"     "0.288"     "0.288"     "0.288"     "0.288"     "0.288"    
[43] "0.288"     "0.288"     "0.288"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 10400"
 [1] "Accepted:" "0.3"       "0.319"     "0.3"       "0.319"     "0.3"      
 [7] "0.319"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[13] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[19] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[25] "0.208"     "0.208"     "0.287"     "0.287"     "0.287"     "0.287"    
[31] "0.287"     "0.287"     "0.287"     "0.287"     "0.287"     "0.287"    
[37] "0.287"     "0.287"     "0.287"     "0.287"     "0.287"     "0.287"    
[43] "0.287"     "0.287"     "0.287"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 10600"
 [1] "Accepted:" "0.3"       "0.319"     "0.3"       "0.319"     "0.3"      
 [7] "0.319"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[13] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[19] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[25] "0.208"     "0.208"     "0.287"     "0.287"     "0.287"     "0.287"    
[31] "0.287"     "0.287"     "0.287"     "0.287"     "0.287"     "0.287"    
[37] "0.287"     "0.287"     "0.287"     "0.287"     "0.287"     "0.287"    
[43] "0.287"     "0.287"     "0.287"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 10800"
 [1] "Accepted:" "0.299"     "0.319"     "0.299"     "0.319"     "0.299"    
 [7] "0.319"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[13] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[19] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[25] "0.208"     "0.208"     "0.286"     "0.286"     "0.286"     "0.286"    
[31] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[37] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[43] "0.286"     "0.286"     "0.286"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 11000"
 [1] "Accepted:" "0.299"     "0.32"      "0.299"     "0.32"      "0.299"    
 [7] "0.32"      "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[13] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[19] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[25] "0.208"     "0.208"     "0.286"     "0.286"     "0.286"     "0.286"    
[31] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[37] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[43] "0.286"     "0.286"     "0.286"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 11200"
 [1] "Accepted:" "0.298"     "0.321"     "0.298"     "0.321"     "0.298"    
 [7] "0.321"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[13] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[19] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[25] "0.208"     "0.208"     "0.286"     "0.286"     "0.286"     "0.286"    
[31] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[37] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[43] "0.286"     "0.286"     "0.286"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 11400"
 [1] "Accepted:" "0.298"     "0.32"      "0.298"     "0.32"      "0.298"    
 [7] "0.32"      "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[13] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[19] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[25] "0.208"     "0.208"     "0.286"     "0.286"     "0.286"     "0.286"    
[31] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[37] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[43] "0.286"     "0.286"     "0.286"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 11600"
 [1] "Accepted:" "0.299"     "0.319"     "0.299"     "0.319"     "0.299"    
 [7] "0.319"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[13] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[19] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[25] "0.208"     "0.208"     "0.287"     "0.287"     "0.287"     "0.287"    
[31] "0.287"     "0.287"     "0.287"     "0.287"     "0.287"     "0.287"    
[37] "0.287"     "0.287"     "0.287"     "0.287"     "0.287"     "0.287"    
[43] "0.287"     "0.287"     "0.287"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 11800"
 [1] "Accepted:" "0.3"       "0.319"     "0.3"       "0.319"     "0.3"      
 [7] "0.319"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[13] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[19] "0.208"     "0.208"     "0.208"     "0.208"     "0.208"     "0.208"    
[25] "0.208"     "0.208"     "0.286"     "0.286"     "0.286"     "0.286"    
[31] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[37] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[43] "0.286"     "0.286"     "0.286"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 12000"
 [1] "Accepted:" "0.301"     "0.319"     "0.301"     "0.319"     "0.301"    
 [7] "0.319"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.286"     "0.286"     "0.286"     "0.286"    
[31] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[37] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[43] "0.286"     "0.286"     "0.286"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 12200"
 [1] "Accepted:" "0.3"       "0.319"     "0.3"       "0.319"     "0.3"      
 [7] "0.319"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.286"     "0.286"     "0.286"     "0.286"    
[31] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[37] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[43] "0.286"     "0.286"     "0.286"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 12400"
 [1] "Accepted:" "0.301"     "0.32"      "0.301"     "0.32"      "0.301"    
 [7] "0.32"      "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.286"     "0.286"     "0.286"     "0.286"    
[31] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[37] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[43] "0.286"     "0.286"     "0.286"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 12600"
 [1] "Accepted:" "0.301"     "0.32"      "0.301"     "0.32"      "0.301"    
 [7] "0.32"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.287"     "0.287"     "0.287"     "0.287"    
[31] "0.287"     "0.287"     "0.287"     "0.287"     "0.287"     "0.287"    
[37] "0.287"     "0.287"     "0.287"     "0.287"     "0.287"     "0.287"    
[43] "0.287"     "0.287"     "0.287"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 12800"
 [1] "Accepted:" "0.302"     "0.32"      "0.302"     "0.32"      "0.302"    
 [7] "0.32"      "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.287"     "0.287"     "0.287"     "0.287"    
[31] "0.287"     "0.287"     "0.287"     "0.287"     "0.287"     "0.287"    
[37] "0.287"     "0.287"     "0.287"     "0.287"     "0.287"     "0.287"    
[43] "0.287"     "0.287"     "0.287"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 13000"
 [1] "Accepted:" "0.303"     "0.32"      "0.303"     "0.32"      "0.303"    
 [7] "0.32"      "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 13200"
 [1] "Accepted:" "0.303"     "0.32"      "0.303"     "0.32"      "0.303"    
 [7] "0.32"      "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 13400"
 [1] "Accepted:" "0.304"     "0.32"      "0.304"     "0.32"      "0.304"    
 [7] "0.32"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 13600"
 [1] "Accepted:" "0.304"     "0.321"     "0.304"     "0.321"     "0.304"    
 [7] "0.321"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 13800"
 [1] "Accepted:" "0.304"     "0.321"     "0.304"     "0.321"     "0.304"    
 [7] "0.321"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 14000"
 [1] "Accepted:" "0.303"     "0.321"     "0.303"     "0.321"     "0.303"    
 [7] "0.321"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 14200"
 [1] "Accepted:" "0.304"     "0.32"      "0.304"     "0.32"      "0.304"    
 [7] "0.32"      "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 14400"
 [1] "Accepted:" "0.304"     "0.32"      "0.304"     "0.32"      "0.304"    
 [7] "0.32"      "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.286"     "0.286"     "0.286"     "0.286"    
[31] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[37] "0.286"     "0.286"     "0.286"     "0.286"     "0.286"     "0.286"    
[43] "0.286"     "0.286"     "0.286"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 14600"
 [1] "Accepted:" "0.304"     "0.319"     "0.304"     "0.319"     "0.304"    
 [7] "0.319"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 14800"
 [1] "Accepted:" "0.304"     "0.319"     "0.304"     "0.319"     "0.304"    
 [7] "0.319"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 15000"
 [1] "Accepted:" "0.304"     "0.32"      "0.304"     "0.32"      "0.304"    
 [7] "0.32"      "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[13] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[19] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[25] "0.212"     "0.212"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 15200"
 [1] "Accepted:" "0.304"     "0.319"     "0.304"     "0.319"     "0.304"    
 [7] "0.319"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[13] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[19] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[25] "0.212"     "0.212"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 15400"
 [1] "Accepted:" "0.303"     "0.318"     "0.303"     "0.318"     "0.303"    
 [7] "0.318"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[13] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[19] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[25] "0.212"     "0.212"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 15600"
 [1] "Accepted:" "0.304"     "0.318"     "0.304"     "0.318"     "0.304"    
 [7] "0.318"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[13] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[19] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[25] "0.212"     "0.212"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 15800"
 [1] "Accepted:" "0.304"     "0.318"     "0.304"     "0.318"     "0.304"    
 [7] "0.318"     "0.213"     "0.213"     "0.213"     "0.213"     "0.213"    
[13] "0.213"     "0.213"     "0.213"     "0.213"     "0.213"     "0.213"    
[19] "0.213"     "0.213"     "0.213"     "0.213"     "0.213"     "0.213"    
[25] "0.213"     "0.213"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 16000"
 [1] "Accepted:" "0.304"     "0.318"     "0.304"     "0.318"     "0.304"    
 [7] "0.318"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[13] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[19] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[25] "0.212"     "0.212"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 16200"
 [1] "Accepted:" "0.304"     "0.317"     "0.304"     "0.317"     "0.304"    
 [7] "0.317"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[13] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[19] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[25] "0.212"     "0.212"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 16400"
 [1] "Accepted:" "0.304"     "0.316"     "0.304"     "0.316"     "0.304"    
 [7] "0.316"     "0.213"     "0.213"     "0.213"     "0.213"     "0.213"    
[13] "0.213"     "0.213"     "0.213"     "0.213"     "0.213"     "0.213"    
[19] "0.213"     "0.213"     "0.213"     "0.213"     "0.213"     "0.213"    
[25] "0.213"     "0.213"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 16600"
 [1] "Accepted:" "0.304"     "0.316"     "0.304"     "0.316"     "0.304"    
 [7] "0.316"     "0.213"     "0.213"     "0.213"     "0.213"     "0.213"    
[13] "0.213"     "0.213"     "0.213"     "0.213"     "0.213"     "0.213"    
[19] "0.213"     "0.213"     "0.213"     "0.213"     "0.213"     "0.213"    
[25] "0.213"     "0.213"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 16800"
 [1] "Accepted:" "0.304"     "0.316"     "0.304"     "0.316"     "0.304"    
 [7] "0.316"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[13] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[19] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[25] "0.212"     "0.212"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 17000"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[13] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[19] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[25] "0.212"     "0.212"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 17200"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[13] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[19] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[25] "0.212"     "0.212"     "0.283"     "0.283"     "0.283"     "0.283"    
[31] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[37] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[43] "0.283"     "0.283"     "0.283"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 17400"
 [1] "Accepted:" "0.302"     "0.315"     "0.302"     "0.315"     "0.302"    
 [7] "0.315"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[13] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[19] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[25] "0.212"     "0.212"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 17600"
 [1] "Accepted:" "0.302"     "0.315"     "0.302"     "0.315"     "0.302"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 17800"
 [1] "Accepted:" "0.302"     "0.315"     "0.302"     "0.315"     "0.302"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 18000"
 [1] "Accepted:" "0.302"     "0.316"     "0.302"     "0.316"     "0.302"    
 [7] "0.316"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[13] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[19] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[25] "0.212"     "0.212"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 18200"
 [1] "Accepted:" "0.302"     "0.316"     "0.302"     "0.316"     "0.302"    
 [7] "0.316"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[13] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[19] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[25] "0.212"     "0.212"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.282"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 18400"
 [1] "Accepted:" "0.302"     "0.316"     "0.302"     "0.316"     "0.302"    
 [7] "0.316"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[13] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[19] "0.212"     "0.212"     "0.212"     "0.212"     "0.212"     "0.212"    
[25] "0.212"     "0.212"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 18600"
 [1] "Accepted:" "0.302"     "0.315"     "0.302"     "0.315"     "0.302"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 18800"
 [1] "Accepted:" "0.302"     "0.315"     "0.302"     "0.315"     "0.302"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 19000"
 [1] "Accepted:" "0.302"     "0.316"     "0.302"     "0.316"     "0.302"    
 [7] "0.316"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 19200"
 [1] "Accepted:" "0.302"     "0.316"     "0.302"     "0.316"     "0.302"    
 [7] "0.316"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 19400"
 [1] "Accepted:" "0.303"     "0.316"     "0.303"     "0.316"     "0.303"    
 [7] "0.316"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 19600"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 19800"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 20000"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 20200"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 20400"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 20600"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 20800"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 21000"
 [1] "Accepted:" "0.303"     "0.314"     "0.303"     "0.314"     "0.303"    
 [7] "0.314"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 21200"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 21400"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 21600"
 [1] "Accepted:" "0.302"     "0.315"     "0.302"     "0.315"     "0.302"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 21800"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 22000"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 22200"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 22400"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 22600"
 [1] "Accepted:" "0.303"     "0.314"     "0.303"     "0.314"     "0.303"    
 [7] "0.314"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 22800"
 [1] "Accepted:" "0.302"     "0.314"     "0.302"     "0.314"     "0.302"    
 [7] "0.314"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 23000"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 23200"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 23400"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 23600"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 23800"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 24000"
 [1] "Accepted:" "0.304"     "0.314"     "0.304"     "0.314"     "0.304"    
 [7] "0.314"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 24200"
 [1] "Accepted:" "0.303"     "0.314"     "0.303"     "0.314"     "0.303"    
 [7] "0.314"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 24400"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 24600"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.283"     "0.283"     "0.283"     "0.283"    
[31] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[37] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[43] "0.283"     "0.283"     "0.283"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 24800"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.283"     "0.283"     "0.283"     "0.283"    
[31] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[37] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[43] "0.283"     "0.283"     "0.283"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 25000"
 [1] "Accepted:" "0.303"     "0.314"     "0.303"     "0.314"     "0.303"    
 [7] "0.314"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.283"     "0.283"     "0.283"     "0.283"    
[31] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[37] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[43] "0.283"     "0.283"     "0.283"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 25200"
 [1] "Accepted:" "0.303"     "0.314"     "0.303"     "0.314"     "0.303"    
 [7] "0.314"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.283"     "0.283"     "0.283"     "0.283"    
[31] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[37] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[43] "0.283"     "0.283"     "0.283"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 25400"
 [1] "Accepted:" "0.304"     "0.314"     "0.304"     "0.314"     "0.304"    
 [7] "0.314"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 25600"
 [1] "Accepted:" "0.303"     "0.314"     "0.303"     "0.314"     "0.303"    
 [7] "0.314"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 25800"
 [1] "Accepted:" "0.304"     "0.313"     "0.304"     "0.313"     "0.304"    
 [7] "0.313"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[13] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[19] "0.211"     "0.211"     "0.211"     "0.211"     "0.211"     "0.211"    
[25] "0.211"     "0.211"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 26000"
 [1] "Accepted:" "0.304"     "0.313"     "0.304"     "0.313"     "0.304"    
 [7] "0.313"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.283"     "0.283"     "0.283"     "0.283"    
[31] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[37] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[43] "0.283"     "0.283"     "0.283"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 26200"
 [1] "Accepted:" "0.304"     "0.313"     "0.304"     "0.313"     "0.304"    
 [7] "0.313"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.283"     "0.283"     "0.283"     "0.283"    
[31] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[37] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[43] "0.283"     "0.283"     "0.283"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 26400"
 [1] "Accepted:" "0.304"     "0.313"     "0.304"     "0.313"     "0.304"    
 [7] "0.313"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.283"     "0.283"     "0.283"     "0.283"    
[31] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[37] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[43] "0.283"     "0.283"     "0.283"     "0.281"    
[1] "Number of not PD:" "0"                
[1] "Iteration: 26600"
 [1] "Accepted:" "0.303"     "0.313"     "0.303"     "0.313"     "0.303"    
 [7] "0.313"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.283"     "0.283"     "0.283"     "0.283"    
[31] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[37] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[43] "0.283"     "0.283"     "0.283"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 26800"
 [1] "Accepted:" "0.304"     "0.313"     "0.304"     "0.313"     "0.304"    
 [7] "0.313"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.283"     "0.283"     "0.283"     "0.283"    
[31] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[37] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[43] "0.283"     "0.283"     "0.283"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 27000"
 [1] "Accepted:" "0.303"     "0.313"     "0.303"     "0.313"     "0.303"    
 [7] "0.313"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.283"     "0.283"     "0.283"     "0.283"    
[31] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[37] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[43] "0.283"     "0.283"     "0.283"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 27200"
 [1] "Accepted:" "0.303"     "0.313"     "0.303"     "0.313"     "0.303"    
 [7] "0.313"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.283"     "0.283"     "0.283"     "0.283"    
[31] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[37] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[43] "0.283"     "0.283"     "0.283"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 27400"
 [1] "Accepted:" "0.304"     "0.313"     "0.304"     "0.313"     "0.304"    
 [7] "0.313"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.283"     "0.283"     "0.283"     "0.283"    
[31] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[37] "0.283"     "0.283"     "0.283"     "0.283"     "0.283"     "0.283"    
[43] "0.283"     "0.283"     "0.283"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 27600"
 [1] "Accepted:" "0.304"     "0.313"     "0.304"     "0.313"     "0.304"    
 [7] "0.313"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 27800"
 [1] "Accepted:" "0.304"     "0.313"     "0.304"     "0.313"     "0.304"    
 [7] "0.313"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 28000"
 [1] "Accepted:" "0.304"     "0.313"     "0.304"     "0.313"     "0.304"    
 [7] "0.313"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 28200"
 [1] "Accepted:" "0.304"     "0.314"     "0.304"     "0.314"     "0.304"    
 [7] "0.314"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 28400"
 [1] "Accepted:" "0.304"     "0.313"     "0.304"     "0.313"     "0.304"    
 [7] "0.313"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 28600"
 [1] "Accepted:" "0.303"     "0.313"     "0.303"     "0.313"     "0.303"    
 [7] "0.313"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 28800"
 [1] "Accepted:" "0.303"     "0.314"     "0.303"     "0.314"     "0.303"    
 [7] "0.314"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 29000"
 [1] "Accepted:" "0.303"     "0.313"     "0.303"     "0.313"     "0.303"    
 [7] "0.313"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 29200"
 [1] "Accepted:" "0.304"     "0.313"     "0.304"     "0.313"     "0.304"    
 [7] "0.313"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 29400"
 [1] "Accepted:" "0.303"     "0.313"     "0.303"     "0.313"     "0.303"    
 [7] "0.313"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 29600"
 [1] "Accepted:" "0.303"     "0.313"     "0.303"     "0.313"     "0.303"    
 [7] "0.313"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 29800"
 [1] "Accepted:" "0.303"     "0.314"     "0.303"     "0.314"     "0.303"    
 [7] "0.314"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 30000"
 [1] "Accepted:" "0.303"     "0.314"     "0.303"     "0.314"     "0.303"    
 [7] "0.314"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 30200"
 [1] "Accepted:" "0.303"     "0.314"     "0.303"     "0.314"     "0.303"    
 [7] "0.314"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 30400"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 30600"
 [1] "Accepted:" "0.303"     "0.315"     "0.303"     "0.315"     "0.303"    
 [7] "0.315"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 30800"
 [1] "Accepted:" "0.302"     "0.314"     "0.302"     "0.314"     "0.302"    
 [7] "0.314"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 31000"
 [1] "Accepted:" "0.303"     "0.314"     "0.303"     "0.314"     "0.303"    
 [7] "0.314"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 31200"
 [1] "Accepted:" "0.303"     "0.314"     "0.303"     "0.314"     "0.303"    
 [7] "0.314"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 31400"
 [1] "Accepted:" "0.303"     "0.313"     "0.303"     "0.313"     "0.303"    
 [7] "0.313"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 31600"
 [1] "Accepted:" "0.303"     "0.313"     "0.303"     "0.313"     "0.303"    
 [7] "0.313"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 31800"
 [1] "Accepted:" "0.303"     "0.313"     "0.303"     "0.313"     "0.303"    
 [7] "0.313"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 32000"
 [1] "Accepted:" "0.302"     "0.313"     "0.302"     "0.313"     "0.302"    
 [7] "0.313"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 32200"
 [1] "Accepted:" "0.302"     "0.313"     "0.302"     "0.313"     "0.302"    
 [7] "0.313"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[13] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[19] "0.209"     "0.209"     "0.209"     "0.209"     "0.209"     "0.209"    
[25] "0.209"     "0.209"     "0.284"     "0.284"     "0.284"     "0.284"    
[31] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[37] "0.284"     "0.284"     "0.284"     "0.284"     "0.284"     "0.284"    
[43] "0.284"     "0.284"     "0.284"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 32400"
 [1] "Accepted:" "0.302"     "0.313"     "0.302"     "0.313"     "0.302"    
 [7] "0.313"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.28"     
[1] "Number of not PD:" "0"                
[1] "Iteration: 32600"
 [1] "Accepted:" "0.302"     "0.313"     "0.302"     "0.313"     "0.302"    
 [7] "0.313"     "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[13] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[19] "0.21"      "0.21"      "0.21"      "0.21"      "0.21"      "0.21"     
[25] "0.21"      "0.21"      "0.285"     "0.285"     "0.285"     "0.285"    
[31] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[37] "0.285"     "0.285"     "0.285"     "0.285"     "0.285"     "0.285"    
[43] "0.285"     "0.285"     "0.285"     "0.28"     
[1] "Number of not PD:" "0"                
